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Abstract 

In this work we apply the Matsubara-Fradkin formalism and the Nakanishi's auxiliary field 
method to the quantization of the Podolsky electrodynamics in thermodynamic equilibrium. This 
approach allows us to write consistently the path integral representation for the partition function 
of gauge theories in a simple manner. Furthermore, we find the Dyson-Schwinger-Fradkin equations 
and the Ward-Fradkin-Takahashi identities for the Podolsky theory. We also write the most general 
form for the polarization tensor in thermodynamic equilibrium. 
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I. INTRODUCTION 



There are several ways to take finite-temperature effects into account in a quantum theory. 
For instance, we can start with the path integral representation for the partition function of 
the system This representation is based on an analogy between the partition function 
and the vacuum-to- vacuum transition amplitude and it is called imaginary-time formalism. 
Once the partition function is evaluated, all thermodynamical quantities can be computed 
from it. The downside of this technique appears most significantly when we consider con- 
strained systems. In those cases we usually tread a long path sorting out the true physical 
degrees of freedom of the theory. For path integrals, this program was developed by Sen- 
janovic j2|. Since all gauge theories are constrained, this step cannot be avoided in this 
approach. Instead of this method, we may consider the real-time formalism, in which the 
time parameter of the field operators becomes a complex number 0, 4|. This method is very 
suitable to study small deviations from equilibrium. Also, in this formalism a natural sepa- 
ration between quantities that depend on the temperature and quantities that do not occurs. 
However, for coherence of the method, the Green functions of the system become matrices, 
with each entry associated with a segment of a path in the complex-time plane. The ap- 
proach known as Thermo Field Dynamics conceives time-evolution of thermal operators as 
well 

a a. 

In this formulation, thermal averages of operators are written as expectation 
values of these operators in the so-called thermal vacuum. However, due to an operation 
called tilde operation, the number of the fields must be doubled. As an alternative to these 
methods, we can begin by defining our ensemble in thermodynamic equilibrium through its 
density matrix. This last formulation was the one designed by Matsubara for non-relativistic 
systems [7']. In its original form, Matsubara's approach only works for systems described 
by the canonical ensemble. Soon after, Fradkin extended Matsubara formalism to include 
also the grand-canonical ensemble and, therefore, he developed a suitable technique to deal 
with (relativistic) quantum field theories, where the number of particles is not fixed js, 9|. 
Fradkin also used Schwinger's source fields in order to make the computation of ensemble 
averages easier. Among the interesting features of the "Matsubara-Fradkin" formalism, we 
mention the following. All positivist quantum theories assume that the state of the system 
is described by a vector in the Hilbert space or, more generally, by a density matrix. In other 
words, it is a postulate of positivist quantum theories that the density matrix of the system 
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contains all the information about the system. So, from the physical point of view, it is more 
natural and fundamental to start from the density matrix. Furthermore, in this method we 
can find a path integral representation for the partition function that is based on no analogy. 
Also, it is possible to extend the formalism to include time-dependent Green functions, de- 
scribing, in this way, small deviations from equilibrium in much the same way the real-time 
formalism does j^. As a final remark, we point out that the spacetime's Euclidean structure 
for systems in equilibrium emerges naturally from the formalism. As a consequence, there is 
a mathematical correspondence between expressions for zero-temperature theories and for 
"hot fields". So, the Matsubara-Fradkin approach to thermodynamical quantization seems 
a reasonable candidate when compared with the other, far more common, techniques. How- 
ever, the Matsubara-Fradkin formalism does not deal, by itself, with the question of gauge 
theories: the quantization of a gauge theory is always a challenge. Even if we consider the 
simplest of all gauge theories, viz. the ordinary electrodynamics, we still face some prob- 
lems in its quantization process. We can take two major alternative routes: we can follow 
Dirac's footsteps or Nakanishi's. In the former case, all the constraints are dealt with in the 
classical level, and we have simple rules relating classical Dirac brackets with commutators 
a.d a,.i-co™™utato. fo. .he ,ua,.u,n theory Q. The main drawback of this method 
is that the Lorentz covariance of the theory is explicitly broken. On the other hand, if 
we follow Nakanishi's quantization procedure, we start with a quantum Lagrangian density 



which incorporates a particular gauge choice ab initio UJ, ll2|]. Due to that choice, the 
theory ceases to be gauge invariant. Despite this explicit breaking of the gauge symmetry, 
Nakanishi's method has one main advantage: it is Lorentz covariant. When we consider a 
gauge theory in thermodynamic equilibrium, there is an advantage in a quantization proce- 
dure that maintains Lorentz covariance in all steps. If we perform the Dirac's constraints 
analysis, in some stage we usually pass from a non-covariant gauge choice to a covariant 
one. This is accomplished using, for example, the Faddeev-Popov trick in the path integral 
representation. By comparison, we see that the Nakanishi's quantization method plays a 
simpler role, being Lorentz covariant from the beginning to the end. If we join together 
Nakanishi's method with the Matsubara-Fradkin formalism we are able to write down the 
correct path integral representations for partition functions of gauge theories without the 
need of performing Senjanovic's procedure for constrained theories. In the present paper 
we deal with the question of quantizing a generalization of the Maxwell electrodynamics in 
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thermodynamic equilibrium. Such a generahzation is known as Podolsky electrodynamics. 

In 1942, Podolsky proposed a theory for the electromagnetic interaction as an alternative 
to Maxwell's 13|. Podolsky theory treats the action for the electromagnetic interaction as 



a functional depending not only on the electromagnetic field and its first derivatives, but 
also on its second-order derivatives. As a consequence, it is possible to split the Podolsky 
field in a sum of a massless vector field with a massive one. Furthermore, Podolsky theory 
depends on a free parameter. This parameter can be interpreted as the mass of the massive 
vector field. Despite the presence of the massive field, Podolsky theory is gauge, as well 
as Lorentz, invariant. As a matter of fact, Cuzinatto, de Melo, and Pompeia proved in 
2007 that Podolsky theory is the only possible extension (regarding theories with second- 
order derivatives) for Maxwell electromagnetism that maintains these two symmetries in 

n 

;he Abelian case [IJ]. Besides, Podolsky electrodynamics can be experimentally tested. In 



15| we, along with co-workers, studied the free Podolsky field in thermodynamic equilib- 



rium using the imaginary-time formalism. In that work we showed that a correction to 
the Stefan-Boltzmann law is expected if the Podolsky field is the underlying field in the 
electromagnetism. Besides, we have used data from cosmic microwave background radiation 
to set a thermodynamical limit to the free parameter of the theory. In 16|], Cuzinatto et 



al. proposed several other experiments for measuring the Podolsky parameter. Moreover, 



m 



17|, one of us with collaborators has shown that the electron self-energy and the vertex 



function of Podolsky theory at zero temperature are both free from ultraviolet divergences 
in the lowest-order of perturbation theory. This does not happen in Maxwellian QED. Our 
goal in this paper is to study Podolsky quantum electrodynamics, also called Generalized 
QED4 (GQED4), in thermodynamic equilibrium. We start with a quick review about the 
Matsubara-Fradkin formalism. Using widely known results from the ordinary QED4, we 
write down the field equations for the fermionic sector of the GQED4. For the gauge sector, 
instead, we use the more convenient Nakanishi's procedure. Once we have all field equa- 
tions, we find the set of equations satisfied by the thermodynamical generating functional of 
Podoslky theory. From the generating functional we obtain the Dyson-Schwinger-Fradkin 
equations and the Ward-Fradkin-Takahashi identities. 

This paper is organized as follows. In section [TTl we review the Matsubara-Fradkin for- 
malism for a general quantum field theory. In section IIIII we study the GQED4 in ther- 
modynamic equilibrium using the Matsubara-Fradkin formalism. In section IIVI we study 
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ensemble averages of some special operators. In section |V] we find path integral representa- 
tions for the thermodynamical generating functional and for the partition function, where 
we show that the Podolsky free partition function coincides with the one we have found in 
our previous work via imaginary-time technique. The Dyson-Schwinger-Fradkin equations 
for the GQED4 in thermal equilibrium are found in section I VII In section IVIII we seek 
the Ward-Fradkin-Takahashi identities for the theory in thermodynamic equilibrium. We 
reserve section IVIIII for our final remarks. As long as Minkowski spacetime is concern, we 
use the metric signature (+, — , — , — ). 

II. THE MATSUBARA-FRADKIN FORMALISM 

In this section we review the Matsubara-Fradkin formalism. As far as we known, although 
Fradkin extensively used this formalism at finite temperature, he only carefully demonstrated 
its zero-temperature version [l^. For this reason, we present his formalism in thermal 
equilibrium from the beginning. 

The density matrix contains all the information about the physical ensemble of the quan- 
tum system. For a system in thermodynamic equilibrium the grand-canonical ensemble 
density matrix reads^ 



where H is the Hamiltonian and is the set of all Noether charges associated with 

continuous internal symmetries of the theory. Each of these charge operators is conserved. 
We assume implicit sum in the index j in this expression over all Noether charges. /3 is 
the inverse of the thermal energy and {fij} is the set of chemical potentials related with 
the Noether charges. We assume H depends on local fields which we denote by 0j (x). 
Since stationarity is a requirement for thermodynamic equilibrium, these fields are time- 
independent. In this section we also assume that our spacetime is (-D + 1) - dimensional. It 
is important to stress that the Hamiltonian H takes into account also all the interactions 
among the fields. 

^ The density matrix also depends explicitly on all the chemical potentials and implicitly on the physical 
(hyper) volume V. We denote it by p (/3) just for the sake of notation. This comment also holds for the 
partition function ([3]). 




(1) 
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The grand-canonical density matrix satisfies the Bloch equation 

dp 

and its trace is the partition function: 



-(H-fi,N,]p{P), (2) 



Z(/3) = Tr[p 



(3) 



We also define the thermal average of an arbitrary operator F as 



Tr 



ZW) 



(4) 



In order to make it easier for us to compute thermal averages we consider the presence of 
classical external sources Sj associated with the several fields The system is, therefore, 
described by the total Hamiltonian Ht- 



Ht = H + Hs 



(5) 



where Hg is the source Hamiltonian. We write it as Hs = J 1-Ls (x) d^x, with the source 
Hamiltonian density l-ig given by 



•H. (x) 



1 



(x)0,(x)-(-lf^ 0,(x) .,(x) 



(6) 



Here, Pj = 1 if the source Sj commutes with its associated field (pj (x) and Pj = 2 if these 
objects anti-commutes with each other, n is the number of fields of the problem. Due the 
presence of the sources, the density matrix now depends on the total Hamiltonian: 



ps{P)= exp 



-13 (Ht - PjNj 



(7) 



We can check that this density matrix satisfies the Bloch equation with H replaced 
by Ht- Since for the great majority of physical problems the fields do not commute with 
H — fijNj, we have ps 7^ P exp (^—f3Hs^ . In order to write p^ in a more convenient 
form, we try the following Ansdtz: 



Ps{P)=p{P)S{P) 
6 



(8) 



Deriving this Ansdtz with respect to /3 and making use of the Bloch equations, we find: 



d/3 



(9) 



Here, we have defined for any operator F its dependence with the temperature in the 
ensemble without external sources: 



F{T)^rHr)Fp{r). (10) 
Integrating ([9]) and using ([I]), (171), and (E]), we get the quantum Volterra equation: 

(11) 



S{/3) = 1- / drH,{r)S{T). 
Jo 

Using an iteration technique, we can find the solution for this equation: 



The ordering operator for two fields is T 











|exp 


- drHs 


(r) 


}■ 










s is T 


'a{t)B{t) 




A{t)B{t) and^ 



(12) 



T 



A (n) B (ra) =9(n- A (n) B (r^) ± 9 {t2 - n) B {r^) A {n] 



(13) 



if Ti 7^ The minus sign is used when both operators are Grassmannian and the plus sign 
otherwise. We define a generalization of the operator S as: 



5(r, r) = T jexp -j dnHs {n] 
Clearly, S 0) = S We can show that the quantity S (r, r') satisfies 



(14) 



S (r, r') = 1 - J dnH, (n) S (ri, r') . (15) 
Now, we functionally derive this expression with respect to the source Sj (x, t^): 



^^^^'^'^ :e{T-T,)9{T,-T')^,{^,T,)S{T,,T')- f d^H , {t,) ^-^j^^ . (16) 



5Sj (x, T^) 



5s j (x, T^) 



^ The Heaviside step function 9{t) is equal to 1 if r > and equal to zero if r < 0. 

7 



Choosing r' = r^. in f[T5]) and multiplying it by (f>j (x, r^.) (xa,, r') results in 



(x, Ta;) (xj;, r') = S (r, r^) 0^ (x, r^) (r^., r) + / driHs (ti) S" (ti, r^.) 0^ (x, r^.) (r^., r') . 

(17) 

Therefore 



SS{t,t') 
6s j (x, r^) 



rfrii^s (ri) 



The solution for this equation is 



9{t - T^) 9 (r^ - r') S (r, r^) 0^ (x, r^^) 

+ 6{t -Tx)6 (r^ - r') / driHg (ri) 5 (ri, r^) 0,- (x, r^) 5* {t^, t') 

J Tx 



5sj (x, r^) ' 



(18) 



SS (r, rQ 
(x, r^) 



: (r - r,) (r,. - r') S (r, r,.) 0, (x, r,.) ^ (r,, r') . 



(19) 



For any operator F, we define its dependence on the temperature in the presence of the 
external sources as 

F^(r)^p7i(r)Fp,(r). (20) 

This definition is very similar to the time evolution of an operator in the Heisenberg 
picture, with the density matrix taking the place of the time evolution operator. For later 
uses, we also compute the derivative of this definition with respect to r:^ 



dF' (r) 
dr 



HT-f^,N„F^iT) 



(21) 



From the very definition of S (r, r'): 



S{t,t,)S{t,,0) = S{t,0) 



(22) 



A,B =AB- BA. 
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Then: 



5S{t,0) 
5s j (x,r^) 



9iT-r,)S{T,0) 05 (x,r,) 



From equations (jH]) and (123!) with < r^. < /3, we get 



5ps (/3) 
5sj (x, r^) 



Ps {(3) (j)] (x, r^) . 



Now, for an arbitrary operator F'^ (r) it is possible to show: 



(5sj (x, r^) 



Hence, we find the result 



(y, Ty) (x, 



Let us define, now, the thermodynamical generating functional: 



Zgf [{Sk}] = Tr [ps (/?)] 



From equation ([3]): 



(23) 



(24) 



(25) 



(26) 



(27) 



ZGF[{sk = 0}] = Zi(3). (28) 

By taking the traces of equations and (!26|) . evaluating the results for the special case 
of vanishing sources, dividing them by the partition function and using (|4]) and fl27|) we find: 



X, Trr 



6 Zgf [{sk}] 



T 



Z (13) 5s j (x, r, 
1 



s=0 
[{Sfc}] 



(29) 
(30) 



Z {P) 5s/(y, Tj^) 5sj(x, r^) 
From these results we see that we can compute thermal averages from the thermodynam- 
ical generating functional. 

Before ending this section, we call attention to the fact that the r variable in f l2^ and 
( 15U]) are restricted to the interval (0, (3). Outside this interval, the functional derivatives can 
be defined in several manners jsl. In section (]VI|) we will define the derivatives outside this 
interval in a consistent way. 



III. GQED4 IN THERMODYNAMIC EQUILIBRIUM 



In this section we consider the GQED4 in thermal equihbrium. We restrict ourselves to 
the (3 + 1) - dimensional case. 

The Podolsky generalized classical electrodynamics is described by the following La- 
grangian density: 

Cged = + ^MC + ^p, (31) 

where 



= i i^a^^,i^b - rrifij^ipa; (32) 

CmC = - QeA^ {Y)ab '^aA; (33) 

Here, 7'^'s are the Dirac matrices satisfying 

{7^7^}., = 2^?^^<5.,, (35) 

V'a = 4'a{x) and = ipai^) c'-re Grassmannian fields, a and b run from 1 to 4, mj is an 
arbitrary parameter with dimension of energy (the so-called "bare fermion mass"), is 
an arbitrary dimensionless parameter (the "bare electrical charge") that plays the role the 
coupling constant of the theory, = Afj,{x) is the Podolsky gauge field, J-'^i, = d^Au — d^A^, 
and mp is the Podolsky parameter, which is real and it has dimension of energy.^ 

We notice that apart from the last term in the right-hand-side of (IMl) . Cged is the 
usual Lagrangian density for electrodynamics. The last term in flMl) is a contribution from 
second-order derivatives of the electromagnetic field 



A. The fermion fields 



As a first step towards a thermodynamical quantum theory, for now we consider only 
Cd + Cmc- This term is exact the same for both Maxwell and Podolsky theories. As it 



The parameter mp can be interpreted as the mass of the massive sector of Podolsky field 15 j 
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is well-known from the Dirac's constraints analysis of Maxwell electrodynamics, the only 
non- vanishing fundamental Dirac brackets between the fermionic fields only is 19 1 



(36) 



Replacing the Grassmannian functions by Grassmannian operators and using the corre- 
spondence principle we get:^ 



-(/).,^(x-y)i. 



(37) 



As we have stated previously, stationarity is a requirement for thermodynamic equilib- 
rium. This means no quantity can depend on time for a system in equilibrium. Since we are 



dealing with fields in thermodynamic equilibrium, we may evaluate ( 1371) . say, for x° 
and simply rewrite it in a time-independent fashion: 







^V^,(x),^,(y)} = -(7°)^,5(x-y)l. 
The grand-canonical density matrix of the GQED4 with external sources is 



Ps {(3) = exp -f3 (JIt - /i, 



where 



(38) 



(39) 



Ht = Hp + Hu + Hmc + Hg] 
1 



Hd=- - I dh\t{^^)^^ (z) , 9f V, (z) 



ruf 



i^a (z),^a (z) 



Hmc = f (7^)., / d'zA, (z) (z) , (z 



H.. 



dh\j^'iz)A^iz) + - r/,(z),^„(z) +- r]a{z),i;^{z) 



(40) 
(41) 

(42) 

(43) 



Hp is a Hamiltonian that commutes with the fermionic fields. J', r/, and rj are the external 
sources. He in ( l39l) is the chemical potential associated with the conserved charge N, which 
reads: 



{Ib} 



AB + BA. 
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(44) 

Now, we apply a similarity transformation with the density matrix to ( !38l) and we use 
fl20l) to write: 

(x, r) , ^ Jy , r) } = - (7°) 5 (x - y ) T. (45) 

This is the fundamental anti-commutation relation for the fermionic part of (Generalized) 
QED in thermodynamic equilibrium. From this and from (12T]) we get the fermionic field 
equations:^ 



N 



(7°) 



ab 



(fz 



i^a (z),^b (z) 



\ 't^ Jab M 
\ 't^ Jba M 



A' 
A' 



- mfSabj i^l (x, r) = r]a (x, r) 1; 
+ mf5ba I (x, r) = (x, r) T, 



(46) 
(47) 



where we have defined the Euclidean Dirac matrices satisfying 



(48) 



and 



Aq =ijP- 



A, 



(49) 
(50) 
(51) 
(52) 



First of all, we notice that the Euclidean spacetime structure of the quantum theory in 
thermodynamic equilibrium has emerged naturally in this formalism. Besides, we see that 
the field equation for t/^ has opposite mass, charge and chemical potential signs if compared 
with the field equation for its conjugated field -0. 



^ From now on, implicit sums on greek indices are taken over an Euclidean geometry. 



12 



B. The Podolsky field 

Now, we apply the Nakanishi's auxiliary field method for the quantization of the Podolsky 
electromagnetic field. In order to do so, we must write a quantum Lagrangian density for 
the Podolsky field which breaks the gauge symmetry explicitly while preserving the Lorentz 
(or, for Euclidean field theories, the SO (4)) symmetry. For the GQED4 in equilibrium, such 
a Lagrangian density is 



A' 



a 



52 ^ ^^^s ^ E^ 
2 2 M V ' M / 



ab 



(53) 



where F^i, = (9^v4^ — d^A^, a is a non- vanishing real parameter (the covariant gauge param- 
eter), B is the auxiliary field, G A is the gauge choice operator, and is the classical 
source.^ By definition, under a gauge transformation, namely. 



A^^A'^ = A, + d,f 
for any well-behaved operator /, the gauge choice operator must satisfy: 



(54) 



G 



A 



G 



^G 



A 



(55) 



Due to this property, the quantum Lagrangian density is not gauge invariant even if 
ge = and = 0. 

The thermodynamical quantum action associated with (!53ll is 



S 



N 



N, 



(56) 



where we have used the notation 



13 Jo Jv 



(57) 



The Schwinger principle for quantum theories states that under an arbitrary infinitesimal 
variation of the field operators the action (l56ll behaves like 



^ From the Euclideanization we have Jo = —ij^ and Jfe = . 
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Sn S'^ — Sn + SSn 



(58) 



with 



6Sn = I d xd^Vf,, 



(59) 



where is a vectorial operator that depends on both the fields and the variations |20|. 
Invoking this principle we find the field equations associated with the Lagrangian density 
dSHD: 



A 



m 



6G* 



A' 



ab 



B = -G 

a 



A' 



where A 



-dfj.dn. Here we have also defined 



(60) 
(61) 



6G 



B- 



A 


6G* 


A 




-h.. = 





SA, 



SA, 



(62) 

where is any vectorial operator and is a suitable vectorial functional of the operators 
B, Af,, and hf,. 

The field equation for the auxiliary field (16T]) is exactly solvable. By replacing the solution 
for B in ( l60l) . its left-hand-side becomes 



6G* 



s 



A' 



SA? 



B 



A + 1") ^s,,A + d,d^) a; ^ 



m 



5G* 



A' 



« 5A? 



■G 



A' 



(63) 



Just like in the Maxwellian case, the differential operator 5^p/S. + df^di, is not invertible. 



A 



m 



The next step in Nakanishi's procedure is finding an appropriate gauge choice G 
order to write the right-hand-side of (163|) as an invertible, local, differential operator acting 
on the gauge field. Inspired by the result found in the classical regime by one of us with a 
collaborator 2l|, we use a gauge choice based on the so-called generalized Lorenz condition: 
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G 



A 



G 



GL 



A 



A 



m 



2 ■ 1 ) 9,A,. 



With this choice, equation ( 1621) becomes 



(64) 



B 



A 



1 d„h 



A 



m. 



+ 1 d^,Bh^ + d^g^, 



where is unimportant for our purposes. This leads us to the identification 



6G* 



A 



SA,, 



With these results, we find the field equation for the Podolsky field: 



Pfiu Al (x, r) 



V^a (X, t) , V'fe (x, r) 



(x,r)l. 



where 



1 / A 
1-- — +1 

a \mp 



dadu 



(65) 



(66) 



(67) 



(68) 



This differential operator differs greatly from its Maxwellian counterpart. For instance, in 
the Maxwell case, the Feynmann-Stiickelberg gauge a = 1 makes the operator independent 
of the term 8^8^. In the Podolsky case, there is no gauge in which that happens. 

C. The ghost fields 



As we have seen, the quantum Lagrangian density (|53|) with the gauge choice f mj) is 
not invariant under an arbitrary gauge transformation fj5^ . However, £jv with = and 
= is invariant under the following special gauge transformation 



^ ^ ^ '■fi — ^ '■fi I 8fj_A, 



A^ A^ — A 
B ^B' = B 



(69) 
(70) 
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provided the operator A satisfies the following constraint 



(A + ml) AA = 0. (71) 

This is the residual gauge symmetry of Cn- It is not difficult to convince yourself that 
this residual gauge symmetry persists in the case Qe 7^ and J^, provided we add the term 
corresponding to free fermions. This is precisely the case of our physical ensemble. 

Now we implement the constraint (ITTIl in the theory by replacing (153|1 by 



Cc = Cn + k\ (x, r) ( ^ + 1 ) AA (x, r) . 

Tfi o 



(72) 



Here, A (x, r) is a Lagrange multiplier operator and k is a constant parameter whose value 
shall be fixed a fortiori. Now we rewrite the Lagrange multiplier as 



A (x, r) = iC (x, r) v. 



(73) 



where C (x, r) is a Grassmannian field operator and t; is a Grassmannian constant. So, up 
to a total derivative of a vectorial operator, we can write 



Cc = Cn + indaC (x, r) 



d 



'^ + S^,]d,C (x,r). 



m 



(74) 



where C 



fields (12I 



^x, r) = vA (x, r) is another Grassmannian field operator. C and C are called ghost 



In order to find the field equations for the ghost fields, we add the Grassmannian fields 
d (x, r) and d (x, r) as sources for them: 



Cgs = Cc + - d,C^ 



d,C 



(75) 



From this quantum Lagrangian density, we find the ghost field equations: 



in (A + 1 ) AC^ (x, T)=d (x, r) 1; 



( A 



1 AC (x,r) 



rffx,r) 1. 



(76) 
(77) 



We have introduced the ghost fields to take the constraint (!7T!) from the residual gauge 
symmetry into account in the GQED4. However, by introducing these news fields, we have 
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incidentally introduced a new symmetry in the theory. The quantum Lagrangian density 
([71]) is invariant under the following global symmetry: 



C fx, t) C fx, 



r)=e^'^°C(x,r); 



Cfx,r) fx,r) =Cfx,r)e 



-ieo 



(78) 
(79) 



where do is a real, constant, arbitrary number. Associated with this continuous global 
symmetry, we have the following conserved charge: 



Q 



- I dh 





+ 




+ 


P,D 


+ 


D,P 



(80) 



where W = 7f(z,r) and tt = 7r(z,r) are the conjugated momentum operators to C and C 
respectively and P = P (z, r) and P = P (z, r) are the conjugated momentum operators to 
D = D (z, r) and D = D (z,r). The extra fields D and D appear from the Ostrogradsky 
formalism for theories of second-order derivatives and are defined in the classical level before 
the Euclideanization as -D = dC/dt and D = dC/dt [22]. They are treated as independent 
of the original ghost fields. The Noether operator Q is called the ghost charge. 

Since we have a new Noether operator in the problem, we must take it into account in 
the density matrix. Therefore, instead of ( 139|) . the true partition function is I23I 



(81) 



Pgs (/3) = exp -13 [Ht- fieN - fIgQ 



where the total Hamiltonian is modified in order to take into account the free Hamiltonian 
for the ghost fields and the source Hamiltonian for the ghost sources. In this expression, 
fig is the ghost chemical potential, i.e., the chemical potential associated with the ghost 
charge operator. We can also define the dependence of an arbitrary operator F with the 
temperature in the ensemble described by this new density matrix. It reads: 



F^'i/3)^Pg,H(3)Fp,sm. (82) 
However, for a field operator wj which commutes with the ghost fields, we have: 
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(/5)= Pas iP) P9s = p7' % P. (/3)= , (83) 

with given by fl39|) . This is the case of the fermionic fields ipa and t/'^ and of the 

Podolsky field and it explains why we can use the density matrix (1391) for these fields. 
However, for the ghost fields themselves: 



C5^(x,r) =e"'^«C^(x,r); (84) 
(x, r) = e-"^«&' (x, r) . (85) 

^ ^ s 

Here, and in field equations fl76|) and f l77|) as well, C*** and C are obtained from a similarity 
transformation with the density matrix (1811) evaluated with fig = 0. 



IV. THE PERIODICITY OF SOME SPECIAL THERMAL AVERAGES 

From the field equations ( l46l) . ( l47l) . ( l67j) . ( j76l) . and ( ITfl) we can find a set of functional 
equations satisfied by the thermodynamical generating functional of the Podolsky theory 
for the electrodynamics. By multiplying each of these field equations by the density matrix 
f lHT]) . making use of fl24|) and fl26l) . taking the trace of the results and using the definition 
(EZD, we find: 



6Z, 



GF 



IK 



A 



m. 



mi 



+ 1 A,. 



X 6r]f, (x, r^) 

5Zgf 
X 6r]b (x, r^.) 

, SZgf 
5Ju (x, r^,.) 



+ 1 A,. 



5d (x, r^.) 
5d (x, r^) 



6' 



^qe (7j) 



5J^{yi,T^)5ri^ (x, 
SJ^{:s.,T^) 6r]b (x, Ta,.) 
(5?7a (x, r^;) (x, Ta;.) 



Va (x, 



+ J„ (x, r^.) 



: - 2{yi,T^) Zgf, 



Zgf', 



(86) 

(87) 

(88) 
(89) 

(90) 
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where Zgf = Zgf [J,V^Vyd.yd] is the thermodynamical generating functionaL The whole 
physical content of the theory is encrypted in this set of coupled functional equations. It 
is possible to solve this set of equations, obtaining in this way the generating functional. 
However, for our purposes, we shall first rewrite this equations in a more convenient form. 
By properly deriving each of these equations once more, we can replace this set of cou- 
pled functional equations by a set of differential equations interconnecting several thermal 
averages of the fields: 



W)a;> ^i''^ - , [5c (y, ry) (x, r.)] ) = ~ ^qe (7^),, 



(^f ^l'''^ + ^f^>^-] , [^'^ (y' ^y) 5^ ^-)] ) 



111/ 



x) T 



(y, Ty) Ay (x, 



IK 



A 



nil 



1 A,.(T 



^«;( +1 ) A,.(T 
mp 



C (y, Ty) C (x, T^) 
C (y, Ty) C (x, r^) 



X (T 



i'c (y> 'Ty) K (x, r^) -il)}, 
+ 5ac5 (x - y) 5 (r^ - Tj/) ; 

^^^^ W)fea 



(91) 



X T 



(y, Ty) Ay (x, r^) ^/-^ 

+ 5ac5 (X - y) 5 (Ta; - Ty) \ 



X, TV 



(92) 



X T 



^« (y> ^y) (X, T^) Ipb 

+ S^^5 (x - y) 5 (r^ - r^^) ; 
: 5 (x - y) 5 (r^. - r^) ; 

: - 5 (x - y) 5 (r^. - Ty) . 



ly, ^s/J 



X, 



where 



Let us consider ensemble averages of the form {T 
is any of the fields. Using the definitions of both the ordering of two operators ( [T3|) 
thermal average (jl]) we write 



(93) 
(94) 

(95) 

V (x, Tx) 

and the 



^r' (y, Ty) (pr (X, 



Tr 



■-d{Ty- T^) 



Pg{(3)(j)r' iy,Ty)(i)r 



X, 



Tr 



±^(r..-r,) 



Pg{l3)^r{^,Tx)^r' (y,' 



(96) 
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where pg is the partition function ( IHT!) evaluated for vanishing classical sources. 

Using an adaptation of equation ffTOj) and the fact that the momentum operator P is the 
generator of the spatial translations, we can write: 



(pr (x, r^) = p^^ (r^) 0r (x, 0) pg (r^) ; (97) 
4)r (x, 0) = e-'^'-^^r (0,0) e'""-^ . (98) 

Now we can compute the traces of the right-hand-side of fl96l) in any basis. Without 
loss of generality, let us consider the mutual eigenbasis for energy, momentum and Noether 
charges \E, P, A^, A^^^). Using this basis for the traces computation, we find that fl^B]) depends 
neither on the absolute values of the spatial coordinates x and y nor on the absolute values 
of the variables related with the temperature Tx and Ty, but only on both the differences 
X — y and Tx — Ty. It is, then, convenient to define the following quantities: 



T Ay (y, Ty) (x, Tx) j = (x - y, r^. - Ty) ; 

T Ty) V'a (X, Tx) ^ = Sab (x - y, Tx - Ty) ] 

T C (y, Ty) C (x, r^) \ =G'(x-y,r^. -Ty). 



(99) 
(100) 
(101) 



From the definition of the ordering operator, we also find 



T 
T 



T 



Af, (X, Tx) Ay (y, Ty) 

A (y, Ty) (x, Tx) 

C (y, Ty) C (X, Tx) 



Duf, (y - X, Tj^ - Tx) ; 

-Sba (y -yi,Ty- Tx) ; 
-G (y - X, - Tx) . 



(102) 
(103) 
(104) 



Since the electromagnetic field commutes with the ghost fields, using f l96|) for < r < /3, 
we find 
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D^t. (x - y, r) 



T 
Tr 



A, (y, 0) (x, t] 
Pg(/3)l^(x,r)l,(y,0) 



Tr 



ZW) 

Pg X (y, /3) (x. ^) 



D^,^ (x - y, r - /3) . 



(105) 



Similarly, taking into account the Grassmannian nature of the fermionic fields and the 
fact that they also commute with the ghost fields, we have: 



Sab (x - y, r) = -Sab (x - y, r - /3) . 



(106) 



So, we see that the thermal averages of the ordering of fields imply the function D^^, must 
be periodic in the r coordinate while the function Sab must be anti-periodic. 

However, repeating these steps for the quantity fllOip and using expressions ( |8^ and (l85ll 
with vanishing sources, we find 



G(x-y,r)= -e-^'^«G(x-y,r-/3) 
G (y - X, -r) = - e^'-^G (y - x, /3 - r) . 



(107) 
(108) 



Equations (El) and ([95]), ffTOTD and ffTMD show us that both G (x - y, - r^) and 
G (y — X, Ty — Tj.) are the Green functions of the operator m (A/m^ + 1) A. This means 
these two quantities must be equal to each other. Hence, it follows: 



G(x-y,r-/3) = e2^'^«G(x-y,r-/3). 



(109) 



Besides the trivial solution pg = 0, the above equation shows us that any ghost chemical 
potential satisfying 



(110) 
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with integer n, is physically admissible. This means that a non-vanishing ghost chemical 
potential is not a real number and, therefore, it is not an observable thermodynamical 
quantity. 

Equation fll07p shows that the ghost Green function G (x, r) can be either periodic or 
anti-periodic in the interval /3, depending on whether we choose n in the chemical potential 
(11 101) to be odd or even, respectively 



23|. 



V. THE PARTITION FUNCTION 

In the present section we will seek path integral representations for the thermodynamical 
generating functional and for the partition function of the theory. In order to achieve that 
goal, we try the following path integral representation for the generating functional: 



JP JA-P J llg 

/ d'^x {J^^Af, + r]J)a - ipaVa - Cd + dC) 
J0 



(111) 



where VA = Y[ T^^a and VipVip = Y[ '^i'a'^i^a- P stands for integration over all periodic 

o-=0 a=l 

configurations of the Podolsky field A^ (x, 0) = A^ (x, /3) and A — P for all anti-periodic con- 
figurations ipa (x, 0) = —ipa (x, /3) and ip^ (x, 0) = —'>Pa (x, /3). These periodicity conditions 
are chosen in order to satisfy the properties ( llOSp for the Podolsky field and ( ]106p for the 
fermionic fields. The integration over the ghost fields C and C depends on the choice for 
the ghost chemical potential (lllOp through (I107p . 

By replacing (II lip in any equation of the set of functional equations (I86tl90p . we find the 
following 

-Zap. (112) 



^3 IX, T) 0(pj (X, T) 



where 0j stands for any of the fields and St = S^''^"''^^ [A, ijj, ijj, C, is the thermodynamical 
action: 
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d 



St = j rf'x j-A^P^V - tnd.C \^5^, + 1 d.C 

The solution for the functional equation (11121) is 



(113) 



Zgf [A^P,^,C,C] =Zoe 



o(/3,Me,V) 



A,i>,ip,C,C\ 



Zn is some so far unknown constant. 



(114) 



By setting the sources equal to zero in equation fillip , we formally get the partition 
function of the GQED4. However, we still have to fix the value of the ghost chemical 
potential in order to perform the path integration over the ghost fields. 

First of all, we notice that the periodicity conditions imposed over the functional integrals 
do not depend whether we are dealing with interacting or with free fields. Taking this 
property into account, we specialize ourselves in the free case, i.e., with = 0. It is possible 
to compute the free partition function exactly: 



Zp {13) = Zo det 



1 M . 



a \m 



det 
p 



A 



ml 



1 A 



det 

Me 



A 



m. 



1 A 



X det 

A-P 



(115) 



Here, detp stands for the determinant evaluated over periodic funcions, detyi-p for the 
determinant evaluated over anti-periodic functions, and det^^ for the determinant over func- 
tions that can be either periodic or anti-periodic depending on the choice of the ghost chem- 
ical potential. From this expression we see that the operator that the ghost determinant 
is evaluated upon appears also in the periodic determinant that comes from the Podolsky 
field, but does not appear in the fermionic, anti-periodic determinant. For this reason, we 
choose the ghost chemical potential f lllOp to be with an odd n [23|. If that is assumed, we 
can rewrite the free partition function as 



Zp{/3) = Zp{/3)Zn{/3) 



(116) 
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where = det^-p 

and 



mf6ab is the partition function for free fermions 



Zp {(3) = det (A + ml) ^ det (A)"^ (117) 

is the partition function for the free Podolsky field, which describes a gas of non-interacting 
massless vectorial particles with Proca particles with mass mp. This result ^incides with 
the one found by us in our previous work with the imaginary-time method [LS]. In writing 
equation f lll7p we have found the values of the two constants that were still unknown, 
namely, k, = oT^I'^ and Zq = detp {mp)~^. 

Using Nakanishi's method and the Matsubara-Fradkin formalism we have found the cor- 
rect partition function for the free Podolsky field without performing its Dirac constraints 
analysis. 

With these new results, we can rewrite the thermodynamical generating functional for 
the interacting theory fillip as 



Zgf = det {mpy^ / VAVCVC / V^jV^je'^^ 



A~P 



xexp 



/ (fx {J^A^+ Va'^a- -iPaVa- Cd + dC) 
J0 



;ii8) 



From this thermodynamical generating functional we can compute any ensemble averages 
of the fields through, e.g., relations (!29|) and (1301) . The Ward-Fradkin-Takahashi identities 
for the GQED4 also follow from this expression. 



VI. THE DYSON-SCHWINGER-FRADKIN EQUATIONS 

In this section we derive the Dyson-Schwinger-Fradkin equations for the GQED4 in ther- 
modynamic equilibrium. The Dyson-Schwinger-Fradkin equations comprise an infinite num- 
ber of coupled nonlinear equations regarding all the Green functions of the theory 18l.l24l.l25|. 

First of all, we notice that the thermodynamical action flllSp is invariant under ipa —'ipa 
and ip^ — )■ —ipa- Therefore, using the solution flllSp for the thermodynamical generating 
functional, we find 
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(V^a(x,r)) = (J„(x,r))=0. (119) 

However, it is not possible to show that the thermal average of the Podolsky field is zero. 
We can only show it is an odd function of the parameter q^. Nevertheless, from the technique 
employed in section HVl we can show that such a quantity is point-independent: 



A^(x,r)) = (A^). (120) 

Now, we define the quantity: 



W [J, 7], r],d,'d\ = In [Zgf [J, r],r],d,d\} . 



(121) 



W is called the thermodynamical generating functional of the connected Green functions. 
We also define the derivatives of this generating functional as: 



Xa (x, r^J 



SW [j, r], r], d, d\ 



SJ^ (x,r^) 
6W [J, r], f], d, d\ 



P[l(x,y;r,,r, 
Sal (x,y;r^.,ry 



6r]a (x, r^) 
8W [ J, 7], T], d, d\ 



(x, r^.) _ 

5Ju (y,7-j/) ' 
5xa (x, r^) 



(122) 
(123) 
(124) 
(125) 
(126) 



^Vh (y, Ty) ■ 

Writing Zgf = in the functional equations ( l86i) and ( 188|) . properly deriving these two 
equations and using these just defined quantities, we find 
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(x - y) (5 (r^ - Ty) = pj^'"'""'* (x, r^.) (x, y; r^., Tj,) - iq^ (7^ ) 



(x,x;r^,r^) 



SJ (x, r^; 



+ ^9 (x, 



(127) 

"^/^acj^iJ (x,y;r^,ry); 

(128) 



-^A. (x, r^.) 



P, 



(nip, a) 



X, r^.) (x, r^.) 



~ (7^ 4a' ^x, r,,) - Xb (x, r^.) Xa (x, 



(129) 



^a(x,r.)= (7f)„,4--) 



5J (x,r^; 



+ (X, 



- mf6ab \ Xb (x,r^) . 



(130) 



Now we define the polarization Il^l and the mass operators implicitly through the 
following relations: 



h")ab ^SfZZT^ - - / ^'^ng (X, z; r., r.) (z, y; r„ r,) ; (131) 

(7j)., ^^ih'''^'"'^ - - / ^'^^i^i 4? (z, y; r., r,) (132) 

With the aid of these two operators, equations (11271) and (11281) show that the quantities 
(11251) and (I126P are the Green functions of the following complete operators: 



(x, y; r., Ty)] - = 6,^6 (x - y) 5 (r.. - r,) > (y) + nW (x, y; r., r,) ; (133) 

52 (x, y; r,, r,)] = (x - y) 5 (r.. - r,) { (7^)^^ D^^^-'^^^ (y, r,)] - m^4.} + 

-SW(x,y;r.,r,). (134) 

The reason why these operators are called "complete" is because they take into account 
all the interactions of the quantum theory. It is worth to emphasize that the functions (I125p 
and (I126P are the Green functions of the operators (I133p and (11340 regardless the presence 
of the external sources. From the very definitions of the polarization and mass operators we 
can prove the following relations 
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(x, y; r^, Ty) = {qef (tJ)^^ / d^ud\ S^} (x, u; r^, r„) rj;' (u, v, y; r„, r^, r^) 

X 5i'J (v, x; r^) ; (135) 
(x, y; r^, ry) = - (ge)^ (7^)^, / d^d'^v Vf^ (x, u; r^, r„) 5^ (x, v; r^, r^) 

J/3 

^ (v, y, u; Tj,, r„) . (136) 

where rj^^^^-j is the complete vertex function: 



^li(a6) y' ^5 ^a;, Ty, T^) — 



(z, r^) 



(7f)^, 5 (x - y) 5 (z-y) 5 (r..-r,) 5 (r, 
^ '^^S (x,y;r^,r^) 

ge 5^9^ (z,r^) 



+ 



We see that in the absence of interaction, T ^(^ab) oc (7^) , . 



(137) 



Computing the complete Green functions ( 1125^ and ( 1126^ for the case of vanishing sources 
and using ([99]), f lToni) . f lTT9l) . and ffT20l) . we find 



(x, y; r^, Ty) = Sab (x - y, - Ty) = Sab (x - y, r^; - Ty) ; (138) 
(x, y; r^, Tj^) = /^/.z. (x - y, - Tj^) - (^A^^ (^A^^ = V^^ (x - y, - Ty) . (139) 

Therefore, in the absence of the classical sources, the polarization and the mass operators 
also depend only on the differences x — y and t^ — Ty. However, as we have called attention 
at the end of section dlT]), the parameter r of the functional derivatives is restricted to the 
interval (0, /3). Due to the periodicity of the function for consistence we shall define the 
functional derivatives outside the original interval, replacing every Dirac delta distribution 
dependent on the r variable that appears in the equation fll33p by the periodic Dirac comb 
distribution (r). A similar statement is true for the delta distributions appearing in the 
equation for Sab, replacing them by the anti-periodic Dirac comb distribution (r). These 
distributions are defined as 
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Aj(r)^ J] {±ir5{T-nP). (140) 

n=— oo 

The Dirac combs satisfy the same periodicity conditions of their corresponding Green 
functions: 



Aj (r - ^) = ±AJ (r) . (141) 

Now, we can show that the complete Green functions in thermodynamic equihbrium 
satisfy: 



(x - y) At (r^ - Ty) = d 



5,,5(x-z)Aj(r.-r,)P,r ^"^ {z) 

+n^a (x - z, - T^)] (z - y, T;, - Ty) ; (142) 
(x - y) A^ (r. - r,) = ! d'z {SacS (x - z) A^ (r, - r,) { (7,^ I^^'^^^ [(l)" 

-rUfScd}^ - ^ad (X - Z, To; - T^) } Sob (z - y, - Ty) . (143) 



In order to clarify the role of the thermal equilibrium, we shall write these expressions 
in Fourier space. Since the quantities appearing in these expressions are either periodic or 
anti-periodic, the Fourier variable associated with r can only assume discrete values: 



Sab (x, 



Aj(r) 



+00 

n=—oo 

+00 

n=— 00 


f d^k ~ 




(144) 

(145) 


^ +00 
n=— 00 






(146) 


+00 
n=—oo 






(147) 



Here, we have defined both the bosonic and the fermionic Matsubara frequencies respec- 
tively as u!^ = 2n7r/P and u;^ = {2n -|- 1) tt/^, for n e Z. With similar Fourier expansions 
for the polarization and mass operators, we find 
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V^4k,u^) =5,,- (148) 



W)ac - ^^^V + <le {K)) - ^Ac - [K^n)] ^cb (k, 0;^) = ^afe- (149) 

~(m? a] 

In these expressions, fc^" = (cu^, k) , fc^" = (^w.f , k) , and P^ii, ^' is the Fourier transform 
of the Podolsky differential operator ( l68l) : 



+ 1 



mi 



1 

a 



+ 1 



ml 



iBriiBn 
i^fi 1^1/ i 1 



(150) 



where (fc^")^ = fc;?"'^;?" = {^uY + k^ > 0. 

We notice that the use of the Dirac comb distributions allowed us to write equations 
for the complete Green functions of the theory in a local fashion in the Fourier space. 
Furthermore, these equations resemble the corresponding expressions at zero temperature. 
The main differences appear in f ll49p . This equation depends explicitly on both the chemical 
potential and the thermal average of the Podolsky operator i^^iij- 

For completeness, we write the equation for the ghost Green function in the Fourier space: 



+ 1 



ml 



(A:^")^GK,k)=l. 



:i5i) 



Despite the ghost fields be Grassmannian, we have fixed the ghost chemical potential 
OllOP to be of the form ^ig = iuj^, which implies a periodicity condition over the ghost Green 
function (11071) . If G is periodic, its Fourier transform must depend on bosonic Matsubara 
frequencies. 

VII. THE WARD-FRADKIN-TAKAHASHI IDENTITIES 

The Ward-Fradkin-Takahashi are relations among the (complete) Green functions of the 



theory 
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In the present section we shall demonstrate them for the GQED4 in ther- 
modynamic equilibrium. 

Let us start by changing the integration functions in the path integral representation for 
the thermodynamical generating functional (11181) accordingly to: 
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(152) 
(153) 
(154) 
(155) 
(156) 

Here, ( (x, r) is a real, periodic function. The generating functional's integration measure 
VADCVCVipVip is invariant under this U{1) transformation 29|. Therefore, the thermo- 
dynamical generating functional can be written as 



A, 


(x, 


r) 


^;(x 


r) 


= (x, r) - (x, r 




(x, 


r) 


^ ^a(x, 


r) 


= e^..C{x,r)^^ (x,r); 




(x, 


r) 


^ V^l(x, 


r) 


= ^,(x,r)e-*^^«^'-); 


c 


(x, 


r) 


^ C7'(x, 


r) 


= C(x,r); 


c 


(x, 


r) 




r) 


= C(x,r), 



Zgf = det {mpy / PAPCDC / VijjVxIje 



a \ m 



1 Ad,A, + d,J, 



' A-P 

X exp <( - ^ d^y |C ' ^ ^ ^ 

-^^.C^ir^'^^^.C - r), (e^"^^ - 1) - (e-^^^« - l) r,. 



xexp 



/ d'^x ( J^A^+ rijja- ipaVa- Cd + dC) 
J0 



(157) 



Since Zgf in equation f lllSp is originally independent of the function ( (x, r), the following 
relation holds: 



5Zgi 



K (x, 

This condition applied to f ll57p implies: 



0. 



(158) 



C=o 



6T 



1 / A 



5r 



-Xa (x, r^j - 



"Xa (X, Tx 



(159) 



.5Xa(x,r^) ' 5xa(x,r^) 
The function F that appears in this expression is the Legendre transform of the functional 



^ We could have defined F as the Legendre transform of the ghost fields as well. However, since they play 
no role in the Ward-Fradkin-Takahashi identities for the Podolsky theory, we have opted for the definition 
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r [d, x,X;d,d\ =W [J, r],T],d,d\ - / d'^x [ J^?9^ + T]^Xa - XaVa] ■ (160) 
By deriving this definition we can show: 



6^r 



5§u (y, Ty) 5df, (x, T^) 



Sal (x,y;r^.,rj, 



-1 



[PW (x,y;r,,r. 



-1 



(161) 
(162) 



From (11591) an infinite number of relations, known as Ward-Fradkin-Takahashi 26N28|. 
follows. For instance, by properly deriving that relation and setting the external sources to 
vanish, we get the following identities: 



5?^r^(,,) (z. 



y, x; r^, Ty, r^) = [Sab (x - z, r^, - r^)] ^ 5 (x - y) (r^ - r. 



[Sab (y - X, - r^.)] ^ 5 (x - z) A 



Ta: - T^. 



1 / A \^ 

[V^^ (x -y,r^ - Ty)]~^ = — + 1 ^9 J (x - y) A+ (r^ - Ty) 

(y. y Tfip J 



(163) 
(164) 



By writing equation (I163p in Fourier space we get the Ward identity in thermodynamic 
equilibrium 



vf^Kab) (k, ; p, uf) = (k + p, u^^i) - Sj (k, ) 
By Fourier-transforming equation (I164p we find: 



(165) 



a 



(fc^")- 



+ 1 



(166) 



Since the right-hand-side of this equation does not depend on the coupling constant 
Qe, this relation holds also in the free case. This can be explicitly checked by comput- 
ing kj^"'P^u ^' (fc^"") from equation (11501) . From equation (I148p we identify the Fourier 
transform of the polarization operator as II^i/ (k, o;^) = (k, w^) — P^ip, ^' (A;-^") and, 
therefore, we have proved the transversality of the polarization operator \2'\ : 
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fc^"n,4k,c^f) =0. (167) 

It is important to emphasize that the transversahty of Uf^^ holds also at finite temperature, 
despite the presence of the Podolsky mass term. The central issue here is that the Podolsky 
theory is gauge invariant, and this symmetry is not broken by thermal effects. 

Using equations ( p9|) and fll39p . we can show the symmetry of the Podolsky Green 
function: P^^, (x, r) = (— x, — r). This implies the symmetry of the polarization: 
Ilf^i, (k, cj^) = IIj,^ (— k, — cj^). From this and from the transversahty condition, we can 
write the most general form for the polarization operator for the GQED4 in thermodynamic 
equilibrium: 



n^. (k,a;f) =A{k,co^) 



l,Bn l,Bn 



h,Bn h,Bn 



B 



+ 



(168) 



where A (k, o;^) = A (— k, —oUn) B (k, w^) = B (— k, —00^). The reason for this un- 
usual form is that an interacting theory in thermodynamic equilibrium comprises a medium. 
In the present case, this medium is a plasma. So, besides the usual momentum fc^", there 
is also the velocity of the medium available for building tensors up. In the plasma rest 
frame, the velocity is proportional to 6^0, which explains the dependence of the polar- 
ization tensor on those quantities. The polarization tensor has the same form in Maxwell 
and Podolsky theories SOj. The only differences are the functions A (k, w,^) and B (k, w^) 
where, in the case of the GQED4, all the dependence on the Podolsky parameter lies. 



VIII. FINAL REMARKS 



Starting from the very concept of density matrix, we have reviewed the Matsubara- 
Fradkin formalism. As we have seen, there is a natural correspondence between the field 
dependence on the temperature-associated variable, denoted by r, through the similarity 
transformation of the field with the partition function with the time evolution of fields at zero 
temperature. Fundamentally, it is precisely this correspondence that justifies the analogy 
between the partition function with the transition amplitude found in the imaginary-time 

32 



formalism. Besides, in this formalism, the Euclidean structure of the spacetime emerges 
naturally. In other words, we do not need to perform a Wick rotation as is usual in other 
techniques. Once we have finished the review of the method, we applied it to the thermo- 
dynamical quantization of the GQED4. 

We have shown that the fermionic part of the electrodynamical theory can be quantized a 
la Dirac without further complications. Since this part of the theory does not depend on the 
explicit form of the Lagrangian density of the gauge field, the quantization presented here 
applies to both Maxwell and Podolsky theories. However, for the gauge part of the theory, 
we have opted for the Nakanishi's covariant quantization. The reason is that even if we had 
done otherwise, we would have had to pass from a non-covariant gauge choice to a covariant 
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one, pretty much as we have done in our previous work in the imaginary time formalism 
We have also seen how the ghost fields arise in the Nakanishi's method and their relation 
with the residual gauge symmetry of the theory. From the field equations for every field 
we have written the set of functional equations satisfied by the thermodynamical generating 
functional. In order to do that, we have made use of the expressions of thermal averages of 
operators in terms of functional derivatives of the thermodynamical generating functional. 
We have proved the periodicity properties of the ordering of the Podolsky and of the fermionic 
fields in thermodyanmic equilibrium. The periodicity of the ghost fields, however, was 
initially depending on the value of the ghost chemical potential, which we showed it must be 
a pure imaginary number and, therefore, not a thermodynamical observable. Then, we have 
found path integral representations for the thermodynamical generating functional and for 
the partition function. Upon computing the partition function of the free Podolsky theory, 
we were able to choose an appropriate value for the ghost chemical potential. The correct 
value for the ghost chemical potential leads to a well-defined periodicity property for the 
ghost Green function. In computing the partition function for the free Podolsky field, we 
have rederived our previous result using a complete different and independent method |l5 |. 
It is important to stress that, in the present formalism, the path integral representation for 
the partition function was found as a solution to the set of functional equations satisfied 
by the thermodynamical generating functional. In other words, we did not need to call 
upon any kind of analogy to construct the correct partition function of the theory. We 
have also written the Dyson-Schwinger-Fradkin equations for the GQED4. Furthermore 
we have written the Ward-Fradkin-Takahashi identities for the GQED4 in thermodynamic 
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equilibrium. From these identities we have proved the Ward identity and the transversahty 
of the Podolsky polarization tensor. Using this result, we have written the most general 
form for the polarization tensor of the GQED4 in thermodynamic equilibrium. Since there 
is a new four-vector available, namely, the four-velocity of the medium, this general form 
differs greatly from its zero temperature version. 

We call attention to the following property of our results: if we take the limit /3 — )■ 00 and 
yU — 7- (which corresponds to the vanishing of temperature and chemical potential) in any 
relation regarding a Green function in thermal equilibrium, we recover the corresponding 
relation in the zero-temperature Euclidean theory. This also holds in the Fourier space. The 
Matsubara frequencies, being them bosonic or fermionic, are inversely proportional to /3. 
In the limit of vanishing temperature, the Matsubara frequencies become dense in the real 
numbers set: they become continuous variables and the sums over them become integrals. 
Furthermore, the Green functions in thermal equilibrium are either periodic or anti-periodic 
with period f3. As this quantity goes to infinity, the thermal Green functions lose their 
periodicity properties and they become usual Euclidean Green functions. As a corollary, if 
we take these limits in any expression involving a Green function, like the Dyson-Schwinger- 
Fradkin and the Ward-Fradkin-Takahashi identities, we find Euclidean versions of those 
found in |l7|]. Besides, since we have not appealed to perturbative techniques, all our results 
hold in both perturbative and non-perturbative regimes of the GQED4. 

Now that we have set the basis for the thermodynamical quantization of the GQED4, 
we can study plasma effects in the GQED4, one of those, for instance, being the Debye 
screening. In the Maxwell theory, the Debye screening consists of the Coulomb electrostatic 
potential becoming short ranged due to the interaction of the field with the plasma. In the 
GQED4 the situation should be even more interesting, since we do not have a simple Coulomb 
potential, but a electrostatic potential a little bit more complicated. Plasma oscillations and 
the fermionic thermal mass generation should also be interesting problems in the GQED4. 
Besides, the Matsubara- Fradkin thermodynamical quantization can be applied, in principle, 
to even more general interactions, like those of non-Abelian gauge theories. 
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